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reference to the thermal expansion properties of various fiber- 
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THERMAL EXPANSION PROPERTIES 
OF COMPOSITE MATERIALS 


Robert R. Johnson, Murat H. Rural, and George B. Mackey 


SUMMARY 


Thermal expansion data for several composite materials, including generic 
epoxy resins, various graphite, boron, and glass fibers, and unidirectional and 
woven fabric composites in an epoxy matrix, have been compiled into one compre- 
hensive report. 

A discussion of the design, material, environmental, and fabrication properties 
affecting thermal expansion behavior is presented. Test methods and their accu- 
racy are discussed. Analytical approaches to predict laminate coefficients of ther- 
mal expansion (CTE) based on lamination theory and micromechanics are also in- 
cluded . 

For space applications, a near- zero CTE is often highly desirable to maintain 
the thermal dimensional stability of a structure on-orbit. Although this is easily 
achievable with composite materials, it is often at the expense of some structural 
efficiency. A discussion is included of methods of tuning a laminate to obtain a 
near-zero CTE. 

Related references that describe these data more extensively are included. 
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INTRODUCTION 

It is well known that the high stiffness and low coefficient of thermal expansion 

V (CTE) of graphite epoxy, together with its low density, make this material espe- 

k 

dally attractive for space applications. 

^ . Much data currently exist as the result of thermal expansion tests on graphite, 

Kevlar, boron, and glass fibers. Properties of composites, including unidirec- 
tional and woven fibers in epoxy , as well as some heat resins , are presented 
here. This document summarizes thermal expansion data obtained from published 
literature, as well as unpublished data from the LMSC data bank. 

Factors that affect the thermal expansion properties are discussed briefly, 
and a related set of references that more extensively describe the influence of 
these factors is presented. 

Use of commercial products or names of manufacturers in this report does not 
constitute official endorsement of such products or manufacturers, either ex- 
pressed or implied, by the National Aeronautics and Space Administration 

FACTORS AFFECTING THE CTE 




The requirement of thermally stable space structures has led to the selection 
of fiber-resin composite materials to achieve a near-zero CTE. A number of fac- 
tors affect the thermal dimensional stability of a laminate, and there is some ques- 
tion as to how close to zero CTE a laminate may be fabricated. Reference 1 pre- 
sents the results of a rather thorough study of 29 test samples from one batch of 
pseudoisotropic GY70/X-30 material. Extension of the CTE data to a large sample 
indicated a 3 sigma statistical variation of close to 0.10 x 10'®/“F (0.18 x 10'®/K). 
Another analytical study (ref. 2), v^^h includes material elastic constants, layup 
angle, material thickness, and input CTE data, indicates that for practical lami- 

nate designs, it is not feasible to specify a CTE less than ±0.05 x 10 /°F 
_6 

(0.09 X 10 /K). This section contains discussions of the following factors 
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which affect the laminate CTE: fiber and void volumes, layup angle, fabric 
skewness, stacking sequence, thermal cycling, temperature dependence, moisture 
effects, and material viscoelasticity. 

Fiber Volume 

The dependence of CTE on fiber volume is illustrated in fig. 1 for a unidirec- 
tional layup. These curves were calculated based on formulas given in ref. 3. As 
seen in this figure , at approximately 60 percent fiber content , the longitudinal 



Figure 1. - Variation of CTE of a undirectional glass 
fiber laminate with fiber volume. 
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CTE is virtually unaffected by any changes in the laminate fiber content . In the 
case of transverse CTE, the sensitivity is more pronounced. In terms of angle ply 
laminates comprised of several layers, the effect of fiber volume variations on the 
thermal expansion behavior of the laminate may not be negligible. 

Void Volume 

The direct effect of voids on the CTE of composite laminates is small within the 
bounds of practical manufacturing requirements (1.5 percent max. void volume). 
However, the presence of voids can indirectly affect the CTE of a laminate by 
initiating microcracks in the resin. Voids in the resin also tend to increase the po- 
tential moisture content of the laminate. Both the microcrack and moisture effects 
are discussed in separate sections. 

Layup Angle 

One of the main advantages of laminated fiber reinforced composites is that mech- 
anical and thermal response of the composites can be tailored directionally to satisfy 
design requirements. This is accomplished by varying the orientation of each layer 
in a systeniatic manner to reach the desired effect. Figure 2 shows the variation 
o^ CTE unidirectional and ±0 angle laminates with fiber orientation (ref. 4). The 
composite CTE can exceed the thermal expansion properties of single layers at cer- 
tain angles for the ±0 angle laminates. This ph-snomenon is predicted by laminate 
theory and also substantiated by tests. 

The sensitivity of the composite CTE to variations from the intended fiber ori- 
entations can be severe. Although manufacturing tolerances for the layup angles 
are typically ±3®, this practice can lead to serious CTE deviations for dimension- 
ally ci’iticnl structures. For example, a 5® deviation in the layup angle for the 
angle ply laminate shown in fig. 2 at approximately 0 = 45° will cause almost an 
order of magnitude change in the value of the CTE of the laminate. Although 
this is an extreme case , it nevertheless points out the degree of sensitivity of 
certain classes of laminates to the change in the layup angle. 
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0, Fiber Direction 


Figure 2. - Variation of expansional strain with laminate fiber orientation 

(from ref. 4). 

Fabric Skewness 

Woven laminae have become quite popular for use with structural composites. 

The effect of fabric skewness is to render the lamina monoclinic (one plane of sym- 
metry). Thus, the lamina will display two axial coefficients and an dependent 
coefficient of thermal shear strain . Such laminae make it impossible to fabricate a 
symmetric laminate. The behavior of unsymmetric laminates is well known. Woven 
laminae should be handled with great care to avoid the inducement of fabric skew- 
ness (e.g. , avoid pulling a fabric on the bias). 

Stacking Sequence 

In general, if a laminate is symmetric and orthotropic (three mutually perpendic- 
ular planes of elastic symmetry), the coefficients of free thermal s^ain are not a 
function of stacking sequence. The presence of coupling response between bending 
4 





and extension behavior in an unsymmetric laminate will generally increase the co- 
efficients of thermal strain (expansion or contraction), because the effective lami- 
nate stiffnesses have been decreased. However, the laminate coefficients of thermal 
curvature are inversely proportional to laminate thickness and can be controlled by 
increasing the number of laminae within the laminate. For coupled laminates, the 
coefficients of thermal strain usually remain quite large compared to their sym- 
metrical counterparts. 

It is well known that if an unsymmetric laminate is restrained to closed geometry 
or zero curvature (e.g., cylindrical tubes or the unsymmetrical facings used to 
fabricate a symmetric laminate sandwich structure) , the general response will ap- 
pear as if the laminate was symmetric. However, for truly sensitive desigi. ap- 
plications involving CTEs, the reliability of this approach is questionable (ref. 5). 

For example, if a tubular member of (Og/+45),p stacking sequence is fabri- 
cated to satisfy a particular design , the tube will not warp , bend , or bow as the 
closed geometry of the tube satisfies the symmetry condition. However, because 
in-plane symmetry conditions are not met in the wall, the tube will exhibit a ther- 
mal shear strain that will take the form of tubular twisting when exposed to ther- 
mal environments. If the tube is proposed to support an antenna, the deleteri- 
ous effect of thermal shear strain is obvious. Thus, for -nsitive thermal expan- 
sion designs , although unsymmetric laminates can be useo , they must be used 
with great care. 


Thermal Cycling 

Thermal cycling is well known to have a significant effect upon the laminat^^ 
CTEs. As shown in refs. 5,6, and 7, the primary influence of thermal cycling is 
to induce resin microcracking. When microcracking occurs, and as it progresses, 
the laminate becomes partially decoupled , both thermally and mechanically . Resin 
degradation proceeds gradually at first, and then somewhat more rapidly, and 
can be detected by changes in the laminate thermal response. For single material 
laminates, the value of the laminate coefficient of expansion drifts toward the uni- 
directional CTE value of the material. This drift depends on many factors, in- 
cluding materials, rate and number of thermal cycles, temperature extremes, 
mechanical load level, and layup angle. For example, the drift of the CTE 



of a cross-pUed HMS laminate with (O/QOg/O)^ layup cycled between -250°F and 
250®F (116 K and 394 K) ' a?, been observed to stop after nine cycles. Laminates 
with more gradual change layer orientations usually take much longer to 
stabilize . 


Temperature Depender e 

It is a common practice to report the CTE of materials as a single quantity; 
these values are often used by designer, and analysts in the same manner. This 
practice may precipitate significant errors in composite design because of the 
temperature dependence of the thermal expansion behavior of composite materials. 
This temperature dependence is mainly caused by the mechanical and physical 
changes in the resin system. For this reason, the CTE values should be obtained 
from thermal expansion test data for the specific design temperature range. The 
CTE is a calculated value which is the slope of the thermal strain-temperature 
curve between two temperatures. This temperature dependence may be observed 
on many of the curves presented later in this document . 

Moisture Effects 

The dimensional stability <^f composites is highly affected by exposure to com- 
plex hj'grothermal histories. Moisture causes swelling and plasticization of the 
resin system. The swelling phenomenon alters internal stresses, thus causing a 
dimensional change in the laminate. Coupled with plasticization of the resin sys- 
tem, ^ .rmanent dimensional changes of up to 30 percent have been observed in 
laboratory experiments after exposure to complex hygrothermal histories (ref. 8). 
The importance of moisture on the thermal expansion behavior of graphite /epoxy 
composites is demonstrated in a study reported in ref. 9. The data shewn in 
table I, :>xtracted from the reference, demonstrate that dimensional changes 
caused by a temperature change of 118®F (339 K) are 2.09 fic, which is much 
less than the dimensional change due to moisture exposure (38.4 Me ). It is there- 
fore important to account for moisture effects in assessing dimensional changes 
due to thermal effects. 
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TABLE I. - CALCULATED STRAINS IN GY70/339 [(O/90)2lg 
LAMINATE DURING HYGROTHERMAL CYCLING 


Sequential hygrothermal history 

Moisture content, 
percent water 

(PO 

After cooling to 75®F (297 K) from 1930F 
(362 K) stress-free temperature 

0.0 

2.09 

After exposure to 75®F (297 K)/95 percent RH 
1 for one year 

0.65 

38.4 

After heating to 160°F (344 K) for 20 minutes 

0.65 

18.9 

After desorbing at 125‘’F (325 K) for 128 days 

0.0 

7.22 

After cooling to 75®F (297 K) 

0.0 

8.11 


Viscoelasticity 


As stated above, changes in the internal stresses due to moisture and thermal 
environment can result in significant dimensional changes in composite laminates. 
Internal stress levels also change due to the viscoelastic phenomenon called the 
relaxation of the resin system. In the true sense, relaxation of the internal stresses 
is a continuous process even without any mechanical, thermal, or moisture excur- 
sions. This pi*ocess is usually minimized by placing fibers in the direction where 
dimensional stability is required. Graphite fibers are known to have negligible 
viscoelastic behavior. 


PREDICTION OF CTE 

The analytical methods used in the prediction of CTEs are presented in the 
following paragraphs. This review includes formulations based on lamination 
theory and the micromechanical approach. A more complete presentation of analyti- 
cal methods discussed here can be obtained from the references listed later in this 
document . 

Successful prediction of laminate thermal properties by the lamination theory is 
highly dependent on defining accurate mechanical and thermal properties of indi- 
vidual layers. In almost all cases, this is accomplished by testing unidirectional 
laminates. The micromechanical approach, where prediction is based on the 
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properties of constituent materials, is seldom used because of several factors. 
These factors include unavailability of constituent material properties, uncertain- 
ties in fiber data, and the inability of the micromechanical approach to define the 
transverse properties of the unidirectional laminate with hifjfh accuracy. 

CTE by lamination Theory 


Classical lamination theory is based upon the concepts of linear anisotropic 
elasticity. Because of the stress and deformation hypotheses tlat are an insepar- 
able part of classical lamination theory , a more correct name would be classical 
thin lamination theory, or even classical laminated plate theory. Strictly speakinfr. 
this theory is valid only for a solid homog^eneous continuum , which is subjected to \ 
homogeneous boundary conditions. The approach and background for this analytical 
method are well developed in refs. 4. 7. and 10 through 12. 

The general constitutive equation for a laminated plate is given by: 


— T 

N + N 


A • B 
1 


m "■ 

T* 

1 

H 

S 

+ 

/ 

i 

1 

B , D 

ij 

K* 


( 1 ) 


This expression is used to relate the laminate mechanical stress resultants and 
stress couples. N. and M. . to the laminate middle surface strains and curvatures, 
t .♦ and K.*. The laminate thermal stress resultants and stress couples arc de- 
fined for a state of uniform (i.e.. homogenous) temperature change: 

T N _K _K 

N. = AT T] Q., [Z - Z„ ,1 Thermal forces (2) 

1 K=1 B 1 n *'■* 

= ~ ^ ^ij Thermal moments (3) 

‘ ‘ K=l 

where . 

' ^OPERATE ' ^REFERENCE* (4) 


*^REFERENCE usually defined as the stress free temperature for the laminate. 
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The laminate atiffhess matrices are also defined : 


*i) = 

N 

E 

K=1 

ii 

'='k 

- *K-l' 

Extensional stiffnesses 

(5) 


1/2 

N 

E 

K=1 


0 2 

iz; • z^., I 

Coupling stiffnesses 

(6) 

“l) = 

1/3 

N 

E 

K=1 

-.I 

tzj - z4_, I 

Bending stiffnesses 

(7) 


The geometric representation of a laminate is depicted in fig. 3. Both the 
laminate and lamina thicknesses (h and hj^l may be defined through the use of the 
following equation: 


h = 


N 

E 

K = l 


‘K 


N 

E 

K=l 


[Z. 


ll 


2 (Y) 



( 8 ) 


- • K — K 

Note that Q., and a . are expressed in terms of the laminate coordinate sys- 
t J I 

tern. The equations necessary to obtain the transformed lamina stiffnesses ) 


r K 


i] 


and transformed coefficients of free thermal strain ( a ^ ) are given in table II. 
These transformations are linear and homogeneous. The invariant forms of the 
equations shown in table II were first given by Tsai and Pagano (ref. 13) in the 



— K — K 

case of Q.. , and Halpin and Pagano (ref. 14) in the case of a • . Note that the 

equations given in table II are sufficient to describe the transformation relation- 
ships for orthotropic laminae. This is a special, but highly practical, case. 

TABLE II. - TRANSFORMATION EQUATIONS FOR AND a 


Term 


FOR THE K"‘ LAMINA 


Rotation 

invariant 



Cos20j^ 

Sin20j^ 

Cos40j^ 

Sin40jj 

«2 

0 

S 

0 

-^2 

0 

u. 

0 

0 

0 

-S 

0 

0 

0 

-«3 

0 

0 

1/2 U 2 

0 

“3 

0 

1/2 U 2 

0 


W 2 

0 

0 

0 

-'^2 

0 

0 

0 

0 

2 W 2 

0 

0 


1/8 

(3Qj^ + 3Q 

22 ^Qi2 

+ 4Q 

1/2 

(Qii 

■ ®22 

) 


1/8 

(Qll 

* ^22 

- 2Qi2 - 

'Q66 

1/8 

^Qll 

*«22 

+ 6Q^2 - 

'Qee 

1/8 

(Qll 

*«22 

- 2Qi2 " 

4^66 

^11 

/(E^l 

-1^2 

^22> 



Wj = 1/2 ( + a^) 

W 2 = 1/2 (Of^ - 


*^11 ^ ^^li 

^22 ^ *"^22 

^12 ^ *^12 ^22 

*^16 = ^26 = ° 









To determine the laminate coefficients of free thermal strain (expansion, 
contraction, or shear) and free thermal curvature (bending or twisting), the me- 
chanical stress resultants and stress couples must be equated to zero: 



(9) 


If equation (9) is now substituted into equation (1) and both sides of the 
resulting equation are divided by ST , the desired solution is obtained : 


T 

II 

'a ! b" 
— 


1 ' 
1 

1 

T 

M 


1 

Q 

PQ 

1 

11 

1 

1 


where : 

♦ 

a . 

1 



Laminate coefficients of free thermal strain 


( 10 ) 


( 11 ) 


* 


= ~T 7 =r Laminate coefficients of free thermal curvature 
1 AT 

Equation (10) can now be used directly and solved by any algebraic method, 
including m Irix inversion. 


Lamina Micromechanics Formulation for CTE 


There are two basic approaches to the micromechanics of composite materials 
(ref. 7): 

(1) Me<'. inics of materials 

(2) ' ’asticity. 

Tue mechanics of materials approach embodies the usual concept of vastly 
s’nplifying assumptions regarding the hypothesized behavior of the mechanical 
system. The elasticity approach may include (1) bounding (variational) tech- 
ni< ues, (2) exact solution, and (3) approximate solutions. 
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Irrespective of the micromechanical approach used, the usual restrictions on 
the composite material are: 

Lamina 

Reinforcement 

Resin 

Macroscopically homogenous 

Homogeneous 

Homogeneous 

Linearly elastic 

Linearly elastic 

Linearly elastic 

Macroscopically orthotropic 
Initially stress- free 

Isotropic 
Regularly spaced 
Perfectly aligned 

Isotropic 


Although the expressions given in the literature for lamina CTEs vary, the 
results of Schapery (ref. 3) are widely quoted. The equation for the longitudinal 
CTE takes the form; 


a 


1 



^m ^m '"m 
¥■ 


m 


m 


( 12 ) 


where : 

Ep = Fiber modulus 

E = Matrix modulus 

m 

= Fiber CTE 
F 

o = Matrix CTE 

m 

'-'c-.v = Volume fractious 

r m 

This equation is quite practical and is suitable for preliminary design. Similarly, 
the transverse CTE can be approximated by (ref. 3) : 


o 


= (1 + 


!’ ) O' 

f' f 




(1 + 


r ) o V 
m m m 


- + 1’ V' ) 

P F F m m^ 


(13) 


However, care must be exercised in using this equation. It is at best an ap- 
proximation . 

References 15 through 20 contain additional information on micromechanical 
methods and CTEs. 
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THERMAL DIMENSIONAL STABILITY AND TUNING FOR ZERO CTE 

This section discusses the thermal dimensional stability of laminated structures , 
including the tuning of laminates to obtain a near- zero CTE with a minimum effect 
upon structural efficiency. 

Thermal Dimensional Stability 

Factors affecting the CTE of a composite were discussed previously. Other 
factors include chemical and physical changes in the resin due to aging and solar/ 
physical radiation. It is significant to note that even during a CTE test under lab- 
oratory conditions, the test specimen is not free from some of these effects, sug- 
gesting that the measured thermal strains in reality are the sum of strains caused 
by different effects. 

Perhaps the most significant design criterion for composite structures requiring 
extreme dimensional stability is the determhiation of the specific mechanical and 
physical environmental spectrum the structure will go through during its lifetime. 
Once the environmental conditions are known, mechanical, physical, and thermal 
behavior of the composite structure can be determined using analytical and experi- 
mental techniques. 

As the starting point , the designer selects the material and laminate configu- 
rations. Stiffness, strength, and CTE requirements determine the type of fi- 
brous materials used in the structure. By determining the number of layers and 
properly orienting them, one can achieve a near- zero CTE laminate that will satisfy 
both the stiffness and strength requirements. Usually, there is more than one 
laminate configuration that will satisfy the preliminary requirements. For better 
thermal dimensional stability, the most likely laminate to succeed is the one that 
is least sensitive to changes affecting the CTE of the laminate. 

A somewhat similar approach has been suggested and shown in refs. 21 and 22. 
A study was conducted on general (O. /±0) laminates to establish a configuration 
with a zero CTE that also satisfies the required mechanical properties. Results 
are shown in fig. 4. 

Laminates with (O/±0)g, and (Og/±0)g layups satisfy the zero CTE 

requirements for a specific 0 . Also, the sensitivity of laminates to ply orientation 
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X 10' 
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[O4/±0ls 

[Og/teJs 


Fi^re 4. - Sensitivity of various laminates to layer onentation 0. 

decreases with an increase in angle. For example, the slope of the expansion 

curve with the (0„/± B) layup is smaller than for the other two laminates men- 

tioned above at zero CTE; thus, less sensitive to changes in angle 9. However, as 

angle 0 gets ’arger, the shear stiffness of the laminate diminishes. To counteract 

this tendency, a 90® ply is added to the (0,/±76) laminate. The resulting lamin- 

ate, (0-/157.5/90) , provides longitudinal extensional stiffness behavior nearly 
0 s 

similar to the (0-/176) laminate, but its shear stiffness is doubled. Results of 
computations for the two laminates are shown below: 


Material and layups 

«x 

E 

_ * 

°xy 



Msi 

(GPa) 

Msi 

(GPa) 

HMS/934 (03/176)^ 

0.229 

(0.41) 

, 16.6 

(112.8) 

1.15 

(7.8) 

HMS/934 (0-/157. 5/90)„ 
•j s 

0.229 

(0.41) 

14.4 

(97.8) 

2.31 

(15.7) 


Once a particular configuration is chosen, further sensitivity studies may be 
performed to establish the dimensional stability of the laminate as influenced by 
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other factors. A typical CTE sensitivity study is shown in table III. In each case, 
the baseline layer properties are increased 10 percent to obtain their effect on the 
composite CTE. 


TABLE III. - (Og/±57.5/90)g LAMINATE CTE SENSITIVITIES TO LAMINA 
PROPERTIES, BASELINE CTE = 0.41 x 10'®/K (0.229 x 10'®/°F) 



Change in laminate CTE 

Lamina properties 

(x 10'®/oF) 

(X 10'®/K) 

^11 

aiCTEl _ 

(10%)Ejj 

-0.10 

^22 

9C-CTE1 = +0 048 

(10%)E22 

+0.09 

^2 

HCTEI 

(10%)Gj2 

0.002 

"l 

^ .m 4 Q Q2g 

0.05 

«2 

^(CTE) ^ g-2 

(10%)O2 

0.09 


Tuning Process 


The tuning of a composite structure to obtain a near-zero CTE with a minimum 
effect on structural efficiency is not an easy task. Most conventional methods to 
drive or tune the CTE of the laminate to zero levels, such as changing the layups 
with different CTE characteristics , will usually result in a sacrifice of structural 
efficiency. 

Theoretically, one of the most effective methods to achieve tuning without 
sacrificing the structural efficiency is to replace a percentage of the laminate 
fibers with an equally stiff but different CTE material. Results of a study show- 
ing the effect of adding boron, silicon carbide (SiC), FP (Al^Og), aluminum, or 
steel fibers to replace a percentage of graphite fibers in an epoxy laminate are 
shown in figs. 5 through 8. Specifically, laminates of VSB-32 pitch and GY 70 
fibers in an epoxy resin were examined. Results are shown in figs. 5 through 6 
for unidirectional laminates and in figs. 7 and 8 for (90/0/ 90) ^ laminates. 

Material properties used in the analysis are given in table IV. 
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Figure 5. - Effect of replacing GY70 fiber with other selected 
fibers in a unidirectional laminate. 
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Figure 6. - Effect of replacing VSB-32 fiber with other selected 
fibers in a unidirectional layup. 
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Figure 7. - Effect of replacing 0° GY70 fiber with other selected 
fibers in a (90/0/90) layup. 



,, Figure 8. - Effect of replacing 0° VSB-32 fiber with other selected 

; fibers in a (90/0/90) layup. 

i 17 

ti 





TABLE IV. - UNIDIRECTIONAL FIBER /EPOXY PROPERTIES - 60 PERCENT FIBER VOLUME 





For GY70 (fig. 5), approximately 20 to 22 percent of the graphite must be 
replaced with boron. SiC, or FP to achieve a zero GTE. This is accomplished 
with almost no reduction in the specific modulus of elasticity for boron and about 
a 5 percent decrease for SiC. 

Only about 7 or 8 percent of boron or SiC is required to obtain the zero CTE 
in unidirectional VSB-32 graphite epoxy (fig. 6). A very small decrease in spe- 
cific stiffness is obtained as a result of using SiC , and a small increase in specific 
stiffness is achieved with boron fibers. 

Steel and aluminum fibers were introduced for comparison purposes. When 
used to replace the VSB-32 fiber, 6 percent of the steel fibe?* required as com- 
pared to 7 or 8 percent with boron or SiC. The decrease in • 'fie modulus of 
elasticity is approximately 9 percent. 

Similar results were obtained for a (90/0/90) laminate in which each of the 
90° plies is 0.003 in. (0.076 mm) thick and the 0° or longitudinal ply is 0.020 in. 
(0.51 mm). For this portion of the study, a percentage of the 0° fibers was re- 
placed by the hybrid fiber. These results are shown in figs. 7 and 8. In fig. 7, 
note that the replacement of about 12 or 13 percent of the GY70 fiber with SiC or 
boron will produce a zero CTE. The specific modulus of elasticity will decrease 

from 563 x lO^ to 551 x 106 \ (143 to 140 : for the boron or SiC 

lbm/xn.3 \ kg/m>>/ 

replacement . 

For the VSB-32 fiber (fig. 8), the CTE of the (90/0/'90) materiel is very close 
to zero. The addition of any replacement fibers will make the laminate CTE posi- 
tive. It should be pointed out that the VSB-32 is still developmental. Any devel- 
opmental increase in fiber modulus will make this (0/90/0) laminate negative; the 
addition of boron or SiC fibers will make the CTE of the laminate zero and alr<^ 

provide an increased stiffness. 

Eo 

An examination of the given in table IV indicates that an efficient proce- 
dure for tuning a laminate with a negative CTE is to select fibers with a high 
positive Ea/p for partial replacement of the 0° fibers, and vice versa. 

Novel Approaches to Tuning 

Some rather innovative techniques have been developed by researchers to 
obtain near-zero CTE laminates, or structural configurations. A tunable end- 
fitting principle for struts is described in ref. 23. 


19 



Another novel approach to the tuning of composite struts is described in 
ref. 24. In particular, a "dual-a'oha" concept is described in detail and a 
"variable-alpha" method is also briefly presented. The dual-alpha strut method 
is a relatively simple concept, which enables the attainment of a specific expansiv- 
ity from parts whose as-cured CTE variability exceeds the design allowables. 

This method uses a laminated construction that exhibits expansh Ities along a 
strut length that are different from one another on either side of a transition zone 
length. In this manner, the expansion characteristics on either side of a transi- 
tion zone bound the desired end-item expansivity. All final tuning <i.e., CTE 
adjustment) is then accomplished in the as-post-cvired condition, thus accounting 
for most material and process variables (e.g. , fiber volume). Tuning is accom- 
plished using the as-cured material only. Addition?’ parts or trning "spacers," 
which are usually of a different material and therefore different thermal diffusiv- 
ity, are not required. This th*'’;i tends to enhance the thermal stability of the 
structure when subjected to transient thermal environments. 

TEST METHODS 

This section presents a brief discussion of the various test methods and 
techniques used to determine the composite CTE. The determination of CTE re- 
quires an approach that uses both length and temperature measurement tech- 
niques. The quality of the CTE data is directly dependent upon the accuracy 
and resolution potential of the equipment utilized. Since composite materials, as 
a class , may CAhibit both very low and very high CTE values , it is often judicious 
to use combinations of test methods in the pursuit of a cost-conscious data- 
generation test plan. Wolff (ref. 25) provides an excellent and comprehensive 
discussion of composite CTE test techniques. Much of this discussion has been 
extracted from this reference. 

Length .Measurement Techniques 

Table V (formulated fn .i refs. 25 and 26) provides a summary of the general 
metrological groups and dimensional measurement techniques in accordance with 
their principles of operation. The accuracy of each method covers wide limits and 
20 
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TABLE V. - SUMMARY OF TECHNIQUES FOR LENGTH MEASUREMENT 


Measurement 

technique 


Resistance strain gages 

Semiconductor strain 
gages 

Capacitance 
LVDT (dilatomr*ers) 
Electronic gages 
Resistance transducers 
Variable impedence 
Variable reluctance 

(1) Non-interferometric 


Resolution, Range, m 
c or (m) or Ae (%) 


Accuracy, m 
or 

range (%) 


A. Electrical transducers 
-7 T 


(O 

20% 

5 X 10'^ 

(O 

1% 

10-® (€) 


o 

1 

10-10 


lo"^ 

o 

1 

00 


10-5 

o 

1 


3 

1% 


10* 2 

0.01% 


10-5 

3% 


B . Electrical — optical transducers 


Contact , 
no contact, 
or small gap. 


Contact 

Contact 
Small gap 
Smedl gap 
Contact 
Contact 
Small gap 
Contact 


! 

Autocollimator 

3 X 10'5 

3% 

10-5 

No contact 

Fiber optics 

10-5 

10^ 

10 "^ 

Small gap 

Optical levers 

1 

o 

io"i 

0.5% 

No contact 

S' igle beam 
shadowing 

2 X 10‘5 

io"5 

0.1% 

No contact 

Scanning beams 

o 

1 

1 

o 

5 X 10 5 

No contact 

Detector arrays 

10-5 

1 -1 
>10 

10-5 ' 

No contact 

Photoelectric 

microscope 

10'^ 

1 

o 

10-5 

i 

Small gap 

Interferometric 

Michelson 

<10"® 

10-5 

10-5 

No contact 

Fabry — Perot 

<10'® 

10-5 

10-5 

Contact 

Fizeau 

10 -» 

10-5 

O 

1 

00 

Contact 

Holographic 

10 "'^ 

10% 

10 -'^ 

No contact 

Speckle 

10 -' 

~1% 

10 -' 

No contact 

Moire 

lO'"^ 

1 

~1% 

1 

1 

o 

No contact 
1 
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TABLE V. - Concluded 


Measurement 

technique 

Resolution, 

(m) 

Range, m 
or Ae (%) 

Accuracy, m 
or 

range (%) 

Contact , 
no contact, 
or small gap, 

<10-1 . 

(m) 

1 

C. Miscellaneous optical 

1 1 1 1 1 



Diffraction pattern 
analysis 

10-® 

>10*^ 

<10''^ 

No 

contact 


Ellipsometry 

lO-'O 

10'5 


No 

contact 


Spatial filtering 

-10“^ 

>10'^ 

I 

o 

No 

contact 


D. Metrological 



-6 

— 


No contact 

Telemicroscopes 

10 

1 

10 


-5 

-2 

-5 


Dial gages 

10 ^ 

10 

10 ^ 

Contact 

Air gage 

lO"^ 

10-3 

1% 

Small gap 


-5 

-1 

..-4 


Micrometers 

3 X 10 ® 

10 

10 

Contact 

Venier scales 

3 X 10*^ 

1 

2 X lO"^ 

Contact 

Profile projectors 

8 X 10“® 

1 

1 

o 

No contact 

Micrometer slides 

10'® 

I 

o 

o 

1 

Small gap 


is, in general, a function of equipment modifications, accessories, and/or the in- 
genuity of the particular designer or user of the apparatus. Several methods suit- 
able for linear displacement measurement on composite materials are described in 
the following paragraphs. Where possible, a primary consideration for a test tech- 
nique is the elimination of contacts to the test sample or structure (e.g. , see 
ref. 25 and table V). 

Electrical Transducers for Displacement Measurement 

Resistance strain gages are attached to the test sample by adhesive bonding. 
These gages are then used to determine the average strain (AL/L^) over a small 
area (i.e. , gage area). Multiple gages can be employed to monitor displacements 
at multiple sites or in multiple directions , and to provide a temperature compensated 
circuit in a Wheatstone bridge unit. This approach is particularly useful for com- 
posite materials if sufficient care is taken to compare identical gage/adhesive sys- 
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terns on a reference material (e.g. , fused silica, in the same bridge circuit). 
Potential problems with the use of these resistance strain gages include gage vari- 
ability, local stiffening of thin sections, moisture, adhesive post cure effects, elec- 
trical heating, lead wires, and contact resistances. Further comments, including a 
brief discussion of semiconductor strain gages, can be obtained from refs. 25 and 
27. 

Capacitance methods, or capacity cell methods, can also be used to determine 
the composite CTE. The principles of operation involve the use of parallel plate 
condensers. The capacitance changes of these plates can be measured as the fre- 
quency changes of an LC oscillator circuit, where the plate spacing is proportional 
to the change in frequency. Calibration is required using identical sample materials 
on the capacitance cell /plate system. This technique is essentially a differential 
one, as described in ref. 25, requiring knowledge of the cell position (upper capac- 
itance plate) during any run. Possible sources of error for capacity cell methods 
include edge capacitance effects, parallelism, sample sizes and shapes, and the time 
required to stabilize temperatures and record readings. 

Linear Variable Differential Transformers (LVDTs) are used widely as a stan- 
dard contact method for determining relative length changes. In practice, voltage 
is induced in the secondary winding of the transformer by the position of the axial 
core of the transducer. Commercial instruments (ref. 25) usually have a dc out- 
put proportional to the axial core position relative to the electrical center. Many 
dilatometers operate with LVDTs as the sensing element that is attached to a 
pushrod placed against an expanding sample in a furnace. Dilatometers using 
LVDTs are marginally useful for low CTE materials because of contacting errors, 
friction, alignment, and temperature effects on the LVDT itself. However. LVDTs 
have potentially infinite resolution, since the signal output is limited only by the 
ability to measure dc voltages, and the vibration stability of the core within the 
LVDT. Sources of error associated with the use of LVDTs and dilatometers, in 
addition to those mentioned above, include the power supply, recording instru- 
ments, and the standard (e.g. . quartz) used to calibrate the LVDT-dilatometer 
system . 

The majority of the LVDT-dilatometer test data given later in this report were 
obtained with the electrical (LVDT) -mechanical (quartz pushrod) system shown in 
fig. 9. The basic approach is described in ref. 26 (ASTM method of vest 228-71). 
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The apparatus shown in fig. 9 uses a specimen 12.00 ±0.010 in. long and 1-in. 
wide . Prior to measurements , the samples are preconditioned in the apparatus for 
1 24 hours at 210 ilO^^F (372 ±6K) with a dry helium purge. Each specimen it then 
subjected to two individual cycles during testing, since it is recognized that CTE 
tends to change with humidity and thermal history. The largest such changes 
usually occur during the first cycle. 

Other dilatometric CTE methods used to generate the LVDT-dilatometer test 
data include the Netzsch dllatometer , modified Leitz dilatometer (fig. 10), ev.^. th. 
dilatometer used in a thermomechaiucul analyzer (TMA). 




Figvire 10. - Modified Leitz dilatometer. 

The Netzsch dilatometer system basically consists of a fused silica stand and 
pushrod between which the test sample is inserted. The movement of the pushrod 
is transmitted to an iron core at the top of the rod that changes the inductance of 
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a transducer , thereby detuning a balanced ac bridge circuit , the rectified 
output of which is proportional to the specimen displacement. 

The modified Leitz dilatometer system shown in fig. 10 is a mechanical- 
optical type method. This instrument is a typical quartz tube dilatometer in which 
the relative distance between ends is transmitted to the recording system by 
quartz rods that contact each end of the specimen. Relative movement of the 
quartz rods is transmitted to a rotatable prism by a mechanical lever system 
that magnifies the specimen displacement. A light beam projected through the 
prism passes through a lens system that further magnifies the displacement. 
Specimen size is restricted only by the size of the dilatometer. An unmodified 
pushrod type Leitz dilatometer is described in ref. 28. 

The use , description , and explanation of the TMA method is well stated in 
ref. 29. This method is not recommended for measurements of low CTE materials, 
as the basic restrictions on specimen size alone may introduce a rather large per- 
centage of error. Finally, electronic gages are often used for gage block compara- 
tors. These gages represent an extension of LVDT-type sensors with exceptional 
care taken to ensure accuracy and repeatability of readings. 

Electrical-Optical Transducers for Displacement Measurement 

Non- Interferometric Methods . - Auto collimaters use collimated light beams to 
detect small angular displacements of a reflecting surface area as small as 

3 

3 x 10 m in diameter. A sensitivity of 0.01 arc sec corresponds to a strain of 

. g 

—4,8 X 10 m/m. This is achievable only with the use of electronic sensors. 

-5 

Visually, only strains of 10 can be detected (ref. 25). These instruments are 
suitable for monitoring the distortions of large structures for straightness or 
flatness . 

Fiber optics bundles provide a noncontact method to provide optical measure- 

.0 

ments over distances of 10 m. The use of fiber optic methods can produce res- 
olutions comparable to LVDTs (ref. 25). An advantage of this approach is that 
light beams can be diverted to reach otherwise inaccessible locations. 



Interferometric Methods . - The interferometric method has been recommended 
by the ASTM (ref. 28). Its upper temperature of application is determined by 
the optical parts of the interferometer (e.g., 1,000 K for vitreous silica). The . 
various interferometric methods are described in refs. 25, 28, 29, 30, and 31. 

An interferometer recombines the two or more parts of a split light beam into 
a third beam whose intensity varies sinusoidally as the optical path length differ- 
ences, relative to one (reference) beam, change by multiples of the wavelength (A). 
The advent of the laser (laser interferometer) has provided the increased coher- 
ence needed to permit the widespread use of this approach for a variety of dimen- 
sional testing applications. 

There are many types of interferometers , but thermal expansion measurements 
mainly employ the original Michelson interferometer (two-beam) , and the Fabry- 
Perot and Fizeau interferometers (multiple reflections). It is fairly easy to count 
intensity variations , or fringes , as a sample expands or contracts when subjected 

_ g 

to thermal excursions. This implies an immediate resolution potential of —10 m. 
For low expansion materials, however, fringe interpolation techniques are re- 
quired. These include analysis of voltage outputs from photodiodes or photo- 
multip.'iers , polarization and phase modulation techniques, smd analysis of photo- 
graphic plates (e.g. , microdensitometry). Values of A/100 - A/1000 are 
presently the state of the art. 

The Michelson interferometry method provides a measure of the relative dis- 
placement of the two reflecting surfaces. Figure 11 (extracted from ref. 25) 
represents a geometrical arrangement applicable to the study of hollow structures 
(e.g., tubes). Fabry-Perot interferometry methods require contact, but offer 
higher resolution potential than the Michelson approach. Fizeau interferometry 
(usually restricted to small samples) normally yields equal thickness circular 
fringes of very large radii of curvature (Newton’s rings) depending on the con- 
tour of the reflecting surfaces. 

Principal disadvantages of interferometric methods include the exposure of the 
measuring equipment to test environments, requirement of highly specialized, ex- 
pensive equipment and difficulties, such as variation of air sources (e.g. , 
changes in index of refraction), and target perturbations which can obliterate the 
time devoted to careful test setups. 


27 



Glass 
belt jar 


Autocollimator 



Tube test 
sample 


Vacuum pressure gage 
*} Liq N 
Helium 

Heating wires 


Strain gage leads 
Thermocouple leads 
T- sensor leads 
Vacuum pump 

MLI (mylar) 
insulation on Cu 
tube with heater 
wires and cooling 
coil 


Laser /detector 


I 

± 
intensity 
stabilization 


Digital 

bidirectional | 
counter 



Frequency 
stabilization I 


Clipper ^ 


Timer 

printer 


Control data 

computer 

analysis 


Damping"| 


Chart 

recorder 


Figure 11. - Laser interferometer and dimensional stability 

test apparatus. 
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Miscellaneous Op 'cal Methods for Displacement Measurement 

Another important method is the diffraction pattern analysis which has been 
used to determine the CTE of E and S glass fibers. A laser beam illuminates the 
test object directly, resulting in a far-field diffraction pattern which is analyzed. 
When the diffracted light is collected with a lens, the pattern formed at its focal 
length is equivalent to the Fourier transform of the test object. Further comments j 
can be found in ref. 25. ^ 

Other methods, described sufficiently elsewhere, are not covered here. 

(See refs. 25, 26, 28, 29, and 30.) In general, the telemicroscope method is 
directly applicable for CTE determinations of fibers , whereas , for example , dial 
gages (coupled with dilatometers as discussed in ASTM D-696, ref. 32) are suit- 
able for the CTE determinations of resins, or very high CTE laminates. 

THERMAL EXPANSION DATA 



This section presents a summary of data relating to the CTE of fibrous com- 
posite materials. The CTE of reinforceifieiits and some resin systems are included 
in addition to the fiber/epoxy data. The data presented were obtained as a result 
of a literature search, and these data have been supplemented with properties 
from the LMSC data bank. Fibers and resins that are no longer available were 
eliminated during the screening process. 

The properties are presented in two formats. Tabular values of the CTE 
are listed for various temperatures as available from the literature. In addition, 
plots of free thermal strain (i.e. , AL/L) as a function of temperature are also 
presented as they are available. In some cases, the data have been replotted 
in an effort to achieve some consistency in the presentation. In other cases, the 
curves have been reproduced as they occurred in the original presentation. 

Tables VI and VII list CTE values for various neat resin systems and fibers, 
respectively. Fiber materials include Kevlar- 49, boron, graphite, and glass. 

Table VIII is a list of available composite thermal test data for unidirectional 
fibers and woven fabrics in epoxy matrix. Each table also includes test methods 
with which the CTE values have been obtained, and the origin of data (i.e. , 
references). 
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TABLE VI. - COEFFICIENT OF THERMAL EXPANSION - 
EPOXY RESIN SYSTEMS 



Temperature 

CTE 



Resin 

OF 

K 

xlO'^T 

xl0“®/K 

Test method 

Source 

Narmco 2387 

RT 

450 

27 

(48.6) 

LVDT 

Ref. 33 


350 


38 

(68.4) 

Dilatometer 


3M PR 279 

RT 

450 

25 

(45) 

LVDT 

Ref. 33 


350 


70 

(126) 

Dilatometer 


Ferro CE-3305 

68-212 

293-373 

35.8 

(64.4) 

LVDT 

Dilatometer 

Du Pont 


Plots of free thermal strain as a function of temperature for various uni- 
directional and fabric composites are given in figs. 12 through 38. These plots 
are useful in determining the CTE of a lamina for a temperature range, 
that is not included in the tables. The value of the CTE is obtained by deter- 
mining the slope of the curves between two temperatures. 






Temperature, K 



(Union Carbide data) 


t 


Figure 12. - Thermal expansion of Kevlar 49/X904B 
(120 style) fabric. 










TABLE VII. - COEFFICIENT OF THERMAL EXPANSION 
FIBROUS REINFORCEMENTS 
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TABLE VIII. - COEFFICIENT OF THERMAL EXPANSION - UNIDIRECTIONAL TAPE 

AND FABRIC EPOXY SYSTEMS 
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Figure 13. - Thermal expansion of unidirectional Boron /5505 tape. 
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(Ref. 43) 


Figure 14. - Thermal expansion of Thornel 300/Narmco 
5208 tape (0*’ fiber dii^ction). 
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Figure 15. - Thermal expansion of Thornel 300/Narmco 
5208 tape (90*^ fill direction). 



Figure 16. - Thermal expansion of HMF 330C/CE339 fabric 
(T-300 fiber) (0° warp direction). 
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Figure 17. - Thermal expansion of HMF 330C/CE339 fabric 
(T-300 fiber) (90° fill direction). 
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(Ref. 36) 


Figure 18. - Thermal expansion of HMF 330C/934 fabric 
(T-300 fiber) (0° warp direction). 


39 


iiCInk'i 





Temperature. K 


280 290 300 310 



Figure 19. - Thermal expansion of HMF 330C/934 fabric 
(T-300 fiber) (90° fill direction). 
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Figure 20. - Thermal expansion of unidirectional 
HTS-2/3501-5A tape (0° fiber direction). 
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Figure 21. - Thermal expansion of unidirectional HTS-2/3501-5A 
tape (90° fiber direction!. 



Figure 22. - Thermal expansion of unidirectional HMS/3501-5A 
tape (0° fiber direction). 
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Figure 23. - Thermal expansion of unidirectional HMS/3501-5A 
tape (90° fiber direction). 
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Figure 24. - Thermal expansion of HMS/3501 unidirectional 

tape (0® direction). 
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Figure 25. - Thermal expansion of KMS/3501 unidirectional 
tape (90° direction). 



(Ref. 45) 
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Figure 26. - Thermal expansion of HMS/CE339 
unidirectional tape (0° fiber direction). 
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Figfure 27. - Thermal expansion of HMS/CE339 unidirectional 
tape (90® fiber direction). 
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Figure 28. - Thermal expansion of HMS/759 unidirectional 
tape (0® direction). 
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Figure 29. - Thermal expansion of HMS/759 
unidirectional tape (90° direction). 
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Figure 30. - Thermal expansion of HMS/934 unidirectional 
tape (0° ftber direction). 
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Figure 31. - Thermal expansion of HMS/934 undirectional 
tape (90° fiber direction). 
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Figure 32. - Thermal expansion of HMS/E788 fabric. 
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Figure 33. - Thermal expansion of unidirectional Thornel-50 
(Pan)/F^63 tape (0® fiber direction). 
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Figrure 34. - Thermal expansion of GY 70 /epoxy 
unidirectional tape fO® fiber direction). 
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Figure 35. - Thermal expansion of GY70/X904B and GY70/934 
unidirectional tape (90° fiber direction). 
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Figure 36. - Longitudinal (0°) thermal expansion 
of P7SS/934, graphite /epoxy . 
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Figure 37. - Thermal expansion of glass /epoxy tooling material 

(181 cloth laminate). 
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Figure 38. - Thermal expartsion of 104 glass scrim cloth. 
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